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■^ ' Abstract 

In general the zero momentum limit of thermal self-energies calcu- 

O ■ lated in perturbation theory depends on the order in which the time 

JL , and the space components of the momentum are taken to zero. We 

QJ \ show that this is an artifact of the perturbative calculation, and in 

i-^ ' fact the non-analyticity of the one-loop self-energy disappears when 

it is calculated with improved vertices and/or improved propagators 
k> i that incorporate the imaginary part of the self-energy. 

H ■ 

The existing calculations of thermal self-energy functions using the for- 
malism of Thermal Field Theory yield results that are not defined if the 
external momentum 4- vector is zero. The classic example is the photon 
self-energy vr^j, in an electron gas. It is well known [|l| that the result of 
the one-loop calculation of 7r^,y(A;°, k) for a photon with external momentum 
k^ = (fc", k) is such that 

Um 7r^,(0,fc)7^ lim7r^,(A;°,0), (1) 

so that the limit in which all components go to zero is not defined. To be 
more specific, this inequality is obtained only for the real part of the self- 



energy. For the imaginary part, the two hniits coincide and hence it is well 
defined. 

There have been two opposite attitudes to this result. On one side it has 
been argued that the inequality should be expected on physical grounds, since 
the quantity lim,^, g7r^i,(0, k) describes the screening of static electric fields 

in the long wavelength limit, whereas limfco^o'^Mi/(^°) 0) gives the plasma 
oscillation frequencies in the quasi-static limit, and they should be different. 
Others believe that this is a mere restatement of the inequality in words 
rather than in symbols, and it does not explain really why the limit of a 
static homogeneous field should be a pathological one. Indeed, Gribosky and 
Holstein even argued, using an effective Lagrangian approach, that the 
limit of zero momentum is a perfectly well-defined one. We will side with 
this second group and consider the inequality of Eq. (0) as a problem, or a 
puzzle. 

Various attempts have been made to resolve this puzzle p|, ||, §, |^, |], 
which involve either introducing new and ad-hoc Feynman rules for thermal 
field theories, or putting restrictions on the general rules 0. Along another 
line of attempt, it has recently been pointed out by Arnold, Vokos, Bedaque 
and Das (AVBD) [§ that the problem mentioned above occurs only if the self 
energy diagram contains two propagators of the same mass. If the masses 
of the particles in the loop are different, the problem does not exist. In 
fact, the calculations of the neutrino self-energy in a gas of electrons and 
nucleons, which were carried out even before the work of Ref . Q] , show this 
feature explicitly P, 0, ^^. It has been speculated that this property may be 



utilized to introduce a mass-splitting regularization for thermal diagrams [Q 
in cases where problems are known to occur. 

This problem, as well the attempts to resolve it, are based on the re- 
sults of one-loop perturbative calculations. It is natural to ask whether the 
singular behaviour of the self-energy function at zero momentum might be a 
consequence of the approximations and idealizations that are implicitly made 
in the perturbative calculations. In this article we show that this is precisely 
the case. 

We begin our analysis by showing that the problem does not exist in 
some theories even if the virtual lines in the loop carry particles of the same 
mass. The example that we give explicitly in this context is the model of a 
fermion ip interacting with a pseudoscalar (p via a Yukawa-type interaction 



of the form ipj^ifjcp. The 1-loop self-energy of the particle in this model is 
well defined in the zero momentum limit. We mention that something simi- 
lar occurs if the interaction is taken to be a derivative coupling of the form 
?/^7^ 75-05^0. While this result is not too surprising since the low momentum 
limit of this theory and the theory with the pseudoscalar coupling is similar, 
this example illustrates in a particularly simple way that the zero momentum 
singularity in some cases disappears if the self-energy diagram is calculated 
with an improved vertex. We then consider some classes of models in which, 
as we show, the 1-loop calculation of the self-energy yields a result that is 
defined at zero momentum when it is carried out using improved propagators 
for the particles that appear in the internal lines of the loop diagrams. Some 
of these propagators have an absorptive part which, as observed by AVBD[|], 
is well-defined at zero momentum even if they are calculated to one-loop and 
the particles in the internal lines of the loop have the same mass. As we will 
see, the absorptive part of such diagrams in turn governs the zero-momentum 
limit of those self-energy diagrams in which the internal lines have the same 
mass. The conclusion that emerges is that, in theories where the problem 
exists, the singular nature of the 1-loop self-energies in the zero momentum 
limit signal a breakdown of the perturbative expansion in that momentum 
regime. In those cases, the calculation must be carried out using an im- 
proved vertex and/or improved propagators in which the absorptive effects 
are taken into account. In a theory with gauge invariance, such as Scalar 
Electrodynamics or QED proper, it is necessary to use both the improved 
vertex and the improved charged-particle propagator since any modification 
to either one of these quantities must be accompanied by a modification to 
the other, in order to satisfy the Ward identity and thereby guarantee gauge 
invariance. 

Before proceeding, we recapitulate some results of the canonical approach 



to the thermal propagators, which we will be using throughout |]I2[. In 
this approach, one has to use anti-time-ordered propagators in addition to 
the time-ordered ones, as well as propagators with no time-ordering. These 
propagators can be arranged in the form of a 2 x 2 matrix. For example, 
for any quantum field ^^ where A denotes any Lorentz index carried by the 



field (none for a scalar field, a Dirac index for the fermion), we can write 

iDt,^{x-y) = (t$^(x)$^(?/)), (2) 



tD^i{x-y) ^ {T <^\x)^\y)) , (3) 

iD^,^{x-y) ^ :F(^\y)^\x)) , (4) 

iDi^{x-y) ^ (<l>^(x)$^(y)) , (5) 

where the bar denotes complex conjugation for bosons and Dirac conjugation 
for fermions and JF is defined by 

_ J — 1 for fermions, , . 

1 1 for bosons. 

T and T are the time-ordering and anti-time-ordering operators defined as 

T$^(x)$''(|/) = Q{x,-y,)<^\x)^''{y)+J'Q{y,-x,)^''{y)^\x)i7) 
T$^(x)$''(|/) = Q{y,-x,)^\x)^''{y)+m{x,-y^)^''{y)^\x)iS) 

where 6 is the step function. We now write the momentum space expansion 
of the field in the form 

^\^) = J j^^ E [ax{p>^{p, A)e-^- + bl{p)v^{p, A)e^^-] , (9) 

where u^ and v^ represent different plane wave solutions arranged by the 
index A, and ax{p) and bx{p) are the annihilation operators for particles 
and antiparticles respectively (for a self-adjoint field ax{p) = h\{p)). The 
properties of the thermal bath come in from the expectation values 

Mp)a\,{p')) = i27rf2E5{p-p')Sxy[l + J'f{p,a)], (10) 

{b,{p)bl,{p')) = {27rf2E6{p-p')S,y[l + J'f{p,-a)], (11) 



with 



f(P^'^)= ^Pp-u-a_jr ^ (12) 



where a plays the role of a chemical potential. We have introduced the 
velocity 4- vector u^ of the heat bath, which has components (1,0) in its own 
rest frame. 

For future purpose, we introduce the matrix U: 



U 



where 



I + J^rj{p,a) Q{—p-u)+J^r]{p,a) 



1 + jr^(p^ a) V 0(P ■ «) + ^V{P, «) 1 + ^V{P, «) 



V{P, «) = 0(P ■ u)f{p, a) + e{-p ■ u)f{-p, -a) 



(13) 



(14) 



for both bosons and fermions, provided we use the appropriate form for the 
function /. In cases where specifically the bosonic or the fermionic form has 
to be used, we will use a subscript B or F to indicate that fact. Thus, the 
procedure described above gives [|1^ the 2x2 scalar field thermal propagator 
as 



A(p) = Ub 



Ao 

-A* 



Uf 



where 



and for fermions 



where 



An = 



p2 - M2 + ZO ' 



sip)^uJl^ _iy, 



(15) 



(16) 



(17) 



^n 



1 



1^ — m + iO ' 

It is easily checked that this gives, for example, the propagators 

1 



(18) 



^uip) = i^ + m^) 
1 



p2 — m'i + iO 



+ 2ni6{p -m^)riF{p,a^) 



An(p) 



p2 - Mj + zO 



27nS{p^~Ml)r]B{p,a^) 



(19) 



for a fermion ifj and scalar 0, which are the ones given by Dolan and Jackiw 



12]. The exphcit forms of the other components are also given in the liter- 



ature |]I1|. We now investigate the implications of these Feynman rules on 
self-energy diagrams in various models. 

Model 1 

Consider the pseudoscalar interaction 

Ant = ^/?^75^$ + A*^ (20) 

Notice that the Lagrangian in this case obeys a parity symmetry 

$^_$, ^^70^, (21) 

under which these are the only possible renormalizable interaction terms. 
A cubic ^^ interaction, for example, is not invariant under this symmetry. 
Using the free-field propagator in Eq. (^) for the fermion field we obtain 

Ren(f(fco,^) = 4/2/-^,^^(p,a^)5(p2-my 



(2^)3 

p ■ k 
k^ + 2p-k 



+ (k^-k)\ , (22) 



which has the unique limit 



RenS)(0,0) = ^fJj^siU + fy)^ (23) 

where we have denoted ^Y f^:z = fF{p,^<^ip) the fermion and antifermion 
momentum distributions. In Eq. (^) the vacuum contribution has been 
omitted, as we will always do henceforth whenever we write explicit expres- 
sions for the self-energies. 

In passing, we note that the property of a unique zero-momentum limit 
does not follow if the fermion-boson interaction is scalar rather than pseu- 
doscalar. On the other hand, if the fermion bilinear is either vector or axial 
vector type, with a derivative coupling to the spin-0 boson, the result of the 
1-loop diagram gives Ren'^]^^(0, 0) = independently of how the limit is 
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taken. Of course, such couplings can arise only if the theory at hand is an 
effective one. Nevertheless, it shows the possibility that the zero momentum 
singularity goes away if the self-energy is calculated with an effective vertex, 
that may itself be the result of an improved calculation, instead of the funda- 
mental one. In the case of the vector coupling, the parity transformation of 
Eq. ( pID is not a symmetry, so it is likely that the theory includes a trilinear 
coupling $^, which gives a contribution to the self-energy that contains the 
zero momentum singularity. 

Model 2 

Consider a scalar $ interacting with another (charged) scalar (f) with the 
following interaction Lagrangian 

Li„t = (/#$* + H.c.) + A^|(^|^ + A.j,|<l>|^ (24) 

It can be easily seen that these are the most general renormalizable interac- 
tion terms if the Lagrangian obeys a global U(l) symmetry under which the 
charges of (j) and $ are 1 and 2 respectively. This symmetry automatically 
rules out any cubic self-interaction of any of the scalar fields. This model is 
similar to the following one 

i^int = (MC^L** + H.C.) + Al*!^ (25) 

in which $ interacts with a Weyl fermion. 

Now consider, for example, the one-loop diagram for the $ self-energy in 
a background of (j) particles, depicted in Fig. |l]. The result of calculating that 
diagram using the free-field propagator given above for the field is 



Ren!f(fco,^) = 2f j j^^r,B{p,aMP 

d?p 



2/ 



(27r)32^ 




k"^ -2p-k k^ + 2p-k 



where p^ = {E, p) with E = Jp^ + M| and we have put /^ ^ = /b(p, ±a</>) to 
denote the momentum distributions. The above formula reveals the problem 
to which we alluded in Eq. (|l]). The same behaviour is obtained for the 
photon self-energy in a background of electrons, or charged scalars. The 
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main observation of this paper is that, as we already mentioned, this problem 
vanishes if the diagram is evaluated employing the full propagator of the (f) 
field instead of the free-field propagator, which is what we show now. 

In what follows we will treat in detail only the two scalar model defined 
by Eq. (|2^) , since the analysis and results are similar for the other model of 
Eq. (pSf). As we mentioned earlier, this problem vanishes if the diagram is 
evaluated employing the full propagator of the field instead of the free-field 
propagator. The full propagator, which we denote by A^'''''{p), can be 
written just hke in Eq. ([T5| ) but with Aq replaced by 

1 

2 ^=i¥y' 



Af = . ,: ^uv (27) 



p2 - Ml - n, 



where IIo is the self-energy function for the (/) field. Thus 


and in particular 







(28) 



Aif (P) = ^ ^ . .^) - 2mp(^\p)MP. «<^) , (29) 



with the spectral density p'-'^-' given by 



imn^^) 



7rp(<^) (p) = -^ ^^ . (30) 

As the interactions are turned off, p'-"'^-' approaches an on-shell delta function. 
However, the existence of a non-zero Im IIo smears the delta function that 
is present in Eq. ([19|) . For this reason it is the easy to see that, if Diagram |l| 
is calculated with the propagator A^^ for the (j) field instead of the free 
particle propagator, then the problem of the non-analyticity of Yin {k) ^t 
zero momentum disappears. Notice that the dispersive part of Hq plays 
no role in this argument. It is not difficult to see that retaining only the 
dispersive part of IIq and neglecting its absorptive part does not remove 



the singularity at zero momentum of Yin {k) 



The next step is to calculate IIq and show that in general it has an 
absorptive part. To this end, we recall that the inverse of the full scalar 
propagator is given by 

A''\p)=p^-M^-U, (31) 

where H is a 2 x 2 matrix whose components must be calculated using the 
Feynman rules of the theory. Comparing Eqs. (^) and (|3T|), the following 
relations are obtained IIT^: 



Hn = no + (no-n*)?7Bb,«) (32) 

U22 = -n* + (no-n*)r^B(p,«) (33) 

ni2 = ~ino-n*)eip-u)fBip,a) (34) 

n2i = -(no-n*)e(-p-w)/B(-p,-a). (35) 

From these, it is easily seen that 

ReUoip) = ReUnip) (36) 

ImHop = — . (37) 

2ifB{p,a) 

Therefore, to determine IIo we must calculate Tlil and TI12 , which can be 
done by evaluating the diagrams in Fig. 0. Since the dispersive part Ho is 
not relevant for resolving the zero momentum problem of the <l> self-energy, 
we will not calculate it here. However, for the consistency of our scheme 
it is important to stress that since the internal lines in Fig. |^ correspond 
to particles of different mass, the function Rell^'j'^ (and hence Re IIq ^) does 
not suffer from the zero momentum problem according to the observation of 
AVBD §. 

We now turn the attention to the calculation of the absorptive part of 
IIq . The simplest way to proceed is to calculate 11^2 &^d then use Eq. (P7|). 
Writing 

k = p' + p, (38) 

the application of the Feynman rules to the diagram of Fig. |^ gives 

- mSt^(p) = i^f)i-^f) J ^^A(ftp')^AS?(fc) , (39) 



where the scalar propagators are given by 

aS?(A;) = -2m6{k^-Ml)fB{k,a^)e{k-u), 

A^^\p') = 2m6{p''-Ml)fBhp',~aMP-^) (40) 

Substituting these propagators into Eq. ( PP| ) and using Eq. ( ^7\) we obtain 

lmUi^\p) = i-nf)eip.u)l0^^Sie-Ml)5ip"-Ml) 

xe{p ■ u)t{k ■ u){fB{p', a^) - fsik, «$)) . (41) 

In writing this formula we have used the identity 

fBi-p\-a^)fBik,a^) = fBip,a^)[fBip\a^) - fBik,a^)] , (42) 

which follows from momentum conservation and the fact that the chemical 
potentials satisfy 

a$ = 2a^ , (43) 

as a consequence of charge conservation. Eq. (^) can be written in the form 

ImU'o^^ = -\p-u\T{p), (44) 

where we have defined 

f t d k u T)^ 

^^^^ " 2^7 (27r)32a;(27r)32^'^^^^ 

X { 6^'\p + p' - k)[U{l + U) - Ml + u)] 
+ s('\p -p'- k)[{l + 4)(1 + U) - f^U)] 

+ 6^'\p + p' + k)[Uf^-{l + U){1 + /^)] 

+ 6^'\p -p' + k)[Ul + f^) - /^(l + /^)] } . (45) 

For the sake of brevity, we have put 

/<^,0 = fB{p',±a^), 

/$,! = fBik,±a^), (46) 

to denote the particle and antiparticle distributions, and in addition, 



p'^ = {E',p'), E' = Jp'^ + Ml, (47) 



k" = {uj,k), uj = V A;2 + Ml . (48) 
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The formula for T given in Eq. ( ^Sj) is immediately recognized as the total 
rate for a (p particle of energy p° and momentum p (as seen from the rest 
frame of the medium) with integrations over the phase space weighted by 
the statistical factors appropriate for each process [|15[. Notice that p is the 



amplitude for the decay processes 0$ — »• cj) and (p -^ 0$, for the annihilation 
processes (p(p —>■ ^ and (p(f)^ — *■ 0, as well as for the inverse reactions of all of 
them. For certain specific values of p^ and p some of these processes will be 
kinematically forbidden, but in general F is non-zero. 

In conclusion, we have presented various models that have similar particle 
content but different interactions, which exemplify several situations. In some 
models the one-loop self-energies are well defined at zero momentum while 
in others the zero momentum singularity dissapears if the coupling is taken 
to be an effective vertex with a suitable momentum dependence. In another 
type of model, the singularity disappears if the self-energy is calculated with 
improved propagators that include the absorptive part. Carrying over this 
idea to the case of QED, which would require that we also use an improved 
vertex in order to maintain the Ward identity and guarantee a gauge invariant 
result, it implies that the one-loop photon self-energy calculated with the full 
charged particle propagator instead of the free propagator is defined at zero 
momentum. In particular, the absorptive part of the charged particle self- 
energy, which physically is related to the damping rate of the particle, cannot 
be neglected if the photon self-energy is evaluated at zero momentum. Then, 
the physical picture that emerges is the following. The traditional formulas 
that are given for 



and 



lim^^,(0,A;) (49) 



hm7r^,(fc°,0), (50) 



which are related to well known physical quantities such as the plasma fre- 
quency and Debye radius, are valid in the limiting cases 

A;0 = ; F<|fc|<m (51) 

k = ; F<A;°<m, (52) 

where m stands for the charged particle mass. Since the two limits correspond 
to two different physical situations the results are different. Traditionally F 

11 



is omitted in the above conditions, but then it must be kept in mind that 
the formulas cannot be taken hterally all the way to zero momentum. 
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Figure 1: Self-energy diagram for a scalar $ in a background of (j) scalar particles in 
the model of Eq. (|2^). The type of each vertex used is depicted in a box near the 
corresponding vertex. 
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Figure 2: Self-energy diagram for cf) in the model of Eq. (|24|) . For the calculation of 
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the left and right vertices should be of type 1 and 2, respectively. 



14 



